Periodic instanton method and macroscopic quantum tunneling between two weakly 

linked Bose-Einstein condensates 
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A new method is used to investigate the tunneling between two weakly-linked Bose-Einstein con- 
densates confined in double-well potential traps. The nonlinear interaction between the atoms in 
each well contributes to a finite chemical potential, which, with consideration of periodic instantons, 
leads to a remarkably high tunneling frequency. This result can be used to interpret the newly 
found Macroscopic Quantum Self Trapping (MQST) effect. Also a new kind of first-order crossover 
between different regions is predicted. 
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Following the first observation of Bose-Einstein con- 
densation (BEC) in dilute gases of trapped alkali atoms, 
remarkable progress has been made both theoretically 
and experimentally [0. In particular, interference be- 
tween two freely expanding condensates has been ob- 
served after switching off the double-well potential that 
confines them . By using a thinner barrier between the 
two condensates it should be possible to establish reliably 
a weak link and study quantum tunneling, or the Joseph- 
son effect, for atoms. Aspects of these questions have 
already been studied theoretically in the limit of nonin- 
teracting atoms ||^ and for small-amplitude Josephson 
oscillations 

Here we develop another theoretical method for a sen- 
sitive and precise investigation of the tunneling between 
two condensates. The almost trivially looking problem of 
the tunneling behavior in a double-well potential has at- 
tracted much attention from theorists for decades. For a 
single particle, the solution can be found even in quantum 
mechanics textbooks Q. The advantage of a nonpertur- 
bative method, as presented here, is that it gives not only 
a more accurate description of the tunneling phenomena 
but also a comprehensive physical understanding in the 
context of quantum field theory. The periodic instan- 
ton configurations, which have been shown to be a useful 
tool in several areas of research such as spin tunneling , 
bubble nucleation and gauge field theory enable 
also the investigation of the finite temperature behavior 
of these systems. In the case of the Bose-Einstein system, 
however, we need to evaluate the tunneling frequency for 
a finite chemical potential even at zero temperature, due 
to the nonlinear interaction between the confined atoms. 
Therefore the chemical potential here replaces the posi- 
tion of the excited energy and gives rise to an expected 
higher tunneling frequency. 

A novel nonlinear effect has been predicted to occur 
in the Bose- Josephson Junction(BJJ) |jl^: The self- 
trapping of a BEC population imbalance arises because 
of the interatomic nonlinear interaction in the Bose gas. 
This was considered to be a novel "macroscopic quantum 
self-trapping" (MQST) and was predicted to be observ- 
able under certain experimental conditions. The three 



parameters, i.e. the ground state energy the inter- 
action energy U , and more importantly, the tunneling 
amplitude K, are still undetermined for a specific ge- 
ometry of the trap and have been taken as constants 
in refs. [ p^|jTl| ]. Here we present a rigorous derivation 
of these quantities and find that they actually depend 
on the number of atoms N . This A^-dependence refines 
the conclusions and makes the self-trapping easier to ob- 
serve. 

The macroscopic wave function $ associated with the 
ground state of a dilute Bose gas confined in the potential 
Vext (r) obeys the well-known Gross-Pitaevskii Equation 
(GPE), which can be obtained using a variational proce- 
dure, i.e. ihd^/dt — SE/S^*. The energy functional E 
is defined by 
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where g = ATrh^a/m is the interatomic coupling constant 
with a the s— wave scattering length. The three terms in 
the integral are the kinetic energy of the condensate Ekm , 
the (an)harmonic potential energy Eho, and the mean- 
field interaction energy Emt , respectively. In the simplest 
case of an isotropic harmonic trap Vexti'r) = 
these energies, which in the Thomas-Fermi approxima- 
tion (TFA) assume the simple values 
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can be calculated beyond the TFA [|l],||jl^] as 
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is a correction term due to 



the presence of a boundary layer near the condensate sur- 
face. Here N is the number of atoms and the harmonic 
oscillator length aho = {h/mujo)^/'^ is introduced for sim- 
plicity. Correspondingly the chemical potential in the 
TFA /iTF = {IhNa/ahof^^ , related to the radius of 
the condensate R through ^tf = raoJ^R^ /2, is modified 
beyond the TFA as ^ = ^tf + 3C/2. We note that in 
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the derivation, the wave function is normahzed to N. If 
one uses instead a wave function normahzed to unity, the 
fohowing correspondence should be reahzed 
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and we obtain the ground state energy E^ ^ ^nd the in- 
teraction self energy f7i,2^i,2 for the isolated traps with 
TVi = 7V2 = as in ref/Q' 
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Considering a condensate of = 5000 sodium atoms 
confined in a symmetric spherical trap with frequency 
Wo = lOOil^, we have E'^ = 1.18nK, UN = 1.03nK, quite 
close to the values estimated in @. 

Calculation of the tunneling frequency by means of the 
periodic instanton method: We study the amplitude for 
tunneling between the two condensates confined in the 
wells of an external double-well potential 



Vextix) 



mojQ 



4r 



(6) 



The two minima are located at zbxo on the a;-axis, and 
the harmonic oscillation frequency near these minima 
is ujQ. The barrier height between the two wells Vq ~ 
is assumed to be high enough so that the over- 
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lap between the wave functions relative to the two traps 
occurs only in the classically forbidden region where in- 
teraction can be ignored and one can safely use the WKB 
wave function approximately The tunneling ampli- 
tude K in ref. can be calculated by different methods, 
and we demonstrate in this work the use of the nonper- 
turbative instanton approach. It is easily shown that this 
tunneling amplitude is just the quantity TZ of pT] (up to a 
minus sign), if one observes the orthogonality property of 
the eigenfunctions / dx^*{x)^j{x) = Sij,i,j = 1,2, with 
*I'i.2(a;) the local modes in each well, which are taken as 
the harmonic oscillator single particle ground state wave 
function in ref. |11|. The nonlinear interaction between 
the atoms in the same well will be included, which modi- 
fies only the chemical potential /i to or beyond the TFA. 

Now we turn to the field theory description of the GPE. 
To this end we consider a scalar field in a 1-dimensional 
time plus 1-dimensional space. After a Wick's rotation 
t — —ir the Euclidean-Lagrangian equation of motion 
for a finite chemical potential takes the form 
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Vext(x) = -fJ, 
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The reason why we can handle a nonlinear problem by 
means of a linear equation of motion is that we discuss 
the tunneling behavior in the barrier region where the 
nonlinear interaction is negligible. However, there are 



obvious differences between the BEC tunneling system 
and the usual one-body problem, i.e. the nonlinear in- 
teraction contributes a finite chemical potential, which 
is just the integration constant on the right hand side of 
eq.(|^). The classical turning points on both sides of the 
barrier can be determined by the relation V{xi_2) — M 
as suggested in ref. For a noninteracting system the 
chemical potential approaches the ground state energy 
corresponding to the vacuum instanton case in . 

Solving this Euclidean time classical equation in the 
usual way [ p^ one obtains the periodic instanton so- 
lution in terms of the Jacobian elliptic function Xc = 
2xokb(k) / luq sn (6(fc)r) with the parameters defined as 
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The Euclidean action for this solution in half of the imag- 
inary period T = 2/C(fc)/6(fc) can be obtained through 
S = /5^/2 {lm{dx/dTf -f Vext{x)) = W + nT/2 with 



W = 



2 8Vo 

3 Wo 



(1 + ^)1/2 (f (fc) - M/C(fc)) 



(9) 



where IC{k) and £{k) are complete elliptic integrals of the 
first and second kinds with modulus fc, respectively. The 
frequency of tunneling between the two condensates is 
then given by the energy level splitting of the two lowest 
states, i.e. = AE/h = 2K/h = 2TZ/h and can be 
calculated by means of the path integral method as |l3[ 
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We emphasize here that this formula has been proven to 
be valid for the entire region when the chemical poten- 
tial is below the barrier height. The condition Vq — fJ. 
determines the sphaleron configuration, where a type of 
phase transition may occur. In the TFA this means 
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where Nx — N1+N2 is the total number of atoms in both 
wells together. Therefore for a specific type of trapped 
atoms and a given double-well potential with separation 
Xq (atom number Nt) there exists a critical number of 
atoms Nci (critical separation Xci) determined by the 
above equation, below (above) which the tunneling pro- 
cess will give the main contribution to the tunneling am- 
plitude. However, above this critical number of atoms or 
below this critical separation value another process, i.e. 
the over-barrier activation will dominate (which is defi- 
nitely not "thermal activation" as in spin tunneling since 
the temperature is zero)(cf. Fig. 1). Between these two 
processes there exists a crossover. A more explicit condi- 
tion for this critical number of atoms (separation between 
the two minima) can be derived beyond the TFA: 
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As an example, we consider two weakly-linked conden- 
sates of Nt — 10** sodium atoms, confined in two sym- 
metric spherical traps witli frequency ujq = 100 Hz as 
in ref. ||l^. The critical value for Xci in the TFA is 
Xci — 24.58/im or more accurately beyond the TFA 
Xci — 25.29/ir7i. We note that here the height of the 
potential barrier is Vb = 2.21nK and the ground state 
is located at hu!o/2 = 0.38riK so that there are several 
energy levels beneath the barrier height. This means the 
interaction between the atoms contributes to the chemi- 
cal potential, which effectively raises the classical turning 
points to a remarkably high level. Although the atoms 
remain in the ground state, the interaction energy is so 
strong that the vacuum instanton method can no longer 
be applied. We have to resort to the periodic instanton 
method, as will be shown below. 



N 



Ncl 



Nc2 



on 




en 






•■g 




_c 


3 




ctivi 


ndo; 


nnel 


uani 


iiddi 


< 


o 


Tu 


O 

u 


Tra 




u 
n 




"5. 


J 






c 


SCO 


Sel 






n 


o 
u 




I 




a 


u 
n 




Ov 











FIG 
parameter 



■^cl ^ c2 

1. Different regions with the corresponding critical 



Low-energy limit: We now consider the "low-energy" 
limit, ^ ^ 0. As in the case of a uniform Bose gas, the 
number of atoms in the ground state can be macroscopic, 
i.e., of the order of the total number in one potential 
well, when the chemical potential becomes equal to the 
energy of the lowest state, which, in our 1-dimensional 
case here, is /ic = ^Tlloq. The lower boundary for the 
chemical potential in fact implies that Rc — aho, i-e. the 
radius of the condensate should never be less than the 
harmonic oscillator length ato- We thus have a result 
similar to that in the vacuum instanton case and the 
"low energy" limit here is only meaningful in this sense. 
Expanding eq. ( |l0|) far below the barrier height, i.e., 
around the modulus /c ^ 1, or equivalently evaluating 
the tunneling amplitude in the vacuum instanton method 
p3[ , we obtain for the tunneling frequency 
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with the Euclidean action 
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This result can be compared with that of ref. [|l| where 
the problem of tunneling in the BEC system is consid- 
ered. We reexpress the tunneling frequency of ref. in 
the notation of the harmonic oscillator length a/jo as 
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which, however, gives not only a smaller exponential con- 
tribution 'SVqI'Hloq (there is a 2/3 factor) but also an inac- 
curate prefactor x\liJ^I2I\ = ujoSc/Ti- The source of this 
inaccuracy is the adoption of the too simple harmonic os- 
cillator wave function of a single particle, which obviously 
oversimplifies the Bose-Einstein condensation tunneling 
problem. At least one should use the WKB wave function 
in the tunneling region, and it can be shown that this cor- 
responds to the vacuum instanton result we present here. 
For the agreement between WKB and vacuum instanton 
methods we refer to ref. |lj] . 

Observation of Macroscopic Quantum Self Trapping: 
The periodic instanton result will lead to a rapidly grow- 
ing behavior for the tunneling frequency when the 
chemical potential, i.e. the number of atoms, is increased. 
According to ref. [|lO| , for a fixed value of the initial pop- 
ulation imbalance z(0) and phase difference 0(0), if the 
parameter A exceeds a critical value Aj,, the population 
becomes macroscopically self-trapped with a nonzero av- 
erage population difference A^i — A^2 • There are different 
ways in which this state can be achieved, and all of them 
correspond to the so-termed MQST condition that 
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This requirement is actually that the modulus of the el- 
liptic function, which appears in the population oscilla- 
tion solution, should be larger than 1 so that the elliptic 
function on will be replaced by dn and the oscillation 
period is shortened from 8klclk)/CA to 4/C(l/fc)/CA. 
The parameters UNt and K are taken as constants in 
ref. jl^ . Considering the fact that they are actually TV- 
dependent as in our calculation above, we can refine the 
conclusions of refs. |lO|JTl| ]. To access the region of self- 
trapping, that is, A > Ac, it is better to lower the value 
of K by making a higher barrier height Vq through in- 
creasing the separation xq or the oscillation frequency 
Wo J than to increase the number of atoms as suggested 
in ref. In fact, the quantity UNt here is propor- 

tional to iiTF ^ N"^/^ which means that increasing the 
number of atoms will not increase the interaction en- 
ergy significantly, and at the same time the tunneling 
amplitude will be increased more drastically. Thus, con- 
trary to the result of refs. 0,0 , we find that the MQST 
will occur when the number of atoms is smaller (instead 
of larger) than a critical value A^c2, i-S- we should de- 
crease the number of atoms instead of increasing it (Fig. 
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1). Inserting the values of the interaction energy and 
the tunnehng amphtude into eq.([l^) we can obtain this 
critical number of atoms for a given potential geometry. 
The parameters which can be adjusted are the number 
of atoms iVr, the oscillation frequency luq, and the sep- 
aration distance between the two condensates xq. Fig 
1. shows the three different regions for different num- 
bers of atoms and distances between the potential wells. 
When Xq (Nt) is smaller (larger) than the critical value 
X2c{Nc2), the atoms will oscillate between these two po- 
tential wells. Once we increase the separation above (or 
decrease the number of atoms below) this critical value, 
the MQST will occur, i.e., most of the atoms will tend to 
remain in their appropriate wells, leading to only a small 
oscillation around a fixed population difference. 

We take the initial condition for the population differ- 
ence to be z{0) = 0.4 and the zero-phase case (j){0) — 
as an example. Other cases with, for example, a non- 
zero phase difference give rise to only a different critical 
parameter Ac. For sodium atoms confined in the double- 
well potential with ujq — lOOHz, we show numerically in 
Fig. 2 the critical line between the three different regions 
in the TFA and beyond it. The upper region marks the 
self-trapping region, the lower the over-barrier activa- 
tion. Quantum tunneling occurs only for a small range 
of the parameter. In the experiment the barrier was 
generated by an off-resonance (blue detuned) laser beam. 
To make our results more applicable to experiment we de- 
note on the right vertical axis the corresponding barrier 
height in units of nK. We also find that the tunneling 
will be suppressed when the separation or the number of 
atoms satisfies xq > 28/im or Nt > 12500 (in the TFA 
Xq > 26^m or Nt > 12000). The crossover will occur 
directly between the self-trapping and the over-barrier 
regions, quite similar to the first-order transition in spin 
tunnehng 0. 
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FIG. 2. Critical line for MQST effect. Solid line: re- 
sults beyond the TFA where the parameters take the values 
UNt = 8/7nTF - 4C and fi = htf + 3C/2. Dashed line: the 
TFA results where UNt = 8/7/iTF and jj, — (itf are used in 
the numerical simulation. 

In conclusion we can say, we have shown that the pe- 
riodic instanton method can be used to investigate the 



tunneling problem in BEG systems at zero temperature. 
The tunneling amplitude and the nonlinear interaction 
energy between the atoms can be calculated analytically 
beyond the TFA. The MQST is more easily observed if 
one takes into account the A^-dependence of the tunnel- 
ing amplitude K and the self interaction energy UNt- 
The crossover between the different regions may be of 
the first-order type when the two minima of the poten- 
tial wells are separated sufficiently far or the number of 
confined atoms in the potential well is large enough. 
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